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I. INTRODUCTION 



Presently there is strong interest to study a newly observed mesons and baryons in the context of a hadronic 
molecule interpretation (for overview see e.g. Ref. In the present work we focus on weak production properties 
of scalar Z)*q(2317) and axial Dsi(2460) charm-strange mesons (for a review see e.g. [^). The Z?*q(2317) meson 
was discovered just a few years ago by the BABAR Collaboration at SLAC in the inclusive Dfn'^ invariant mass 
distribution from e+e^ annihilation data Q. The nearby state Z?si(2460) decaying into D^n^ was observed by the 
CLEO Collaboration at CESR . Both of these states have been confirmed by the Belle Collaboration at KEKB [H] . 
In the interpretation of these experiments it was suggested that the D*g(2317) and Dsi(2460) mesons are the P-wave 
charm-strange quark states with spin-parity quantum numbers = 0"*" and — 1+, respectively. In the following 
the Belle and BABAR [7} Collaborations observed the production of £>*o(2317), I?si(2460) and their subsequent 
strong and radiative transitions in the nonleptonic two-body B decays. The most recent data of BABAR and Belle 
on two-body i?-meson decays into -D*q(2317) and Dsi{2460) states can be found in Refs. [1, It is worth noting 
that the existing experimental information on the properties of Z)*q(2317) and _Dsi(2460) mesons leaves quite a 
significant uncertainty in their interpretation as = 0^ and — 1^ states. 

Theoretical analysis oi B ^ Z?^*^L'*q(I?si) decays has been performed in different approaches [ll|-[23| based on 
the factorization hypothesis, which essentially simplifies the calculation of the transition amplitude. The factorizable 
amplitude B D'-*Wla{Dsi) is given by the product of the corresponding form factors (or their combination) 
describing semileptonic B D^*^lDi transitions and the leptonic decay constant /n*^ {fo^i)- The leptonic decay 
constants fw and fosi have been calculated directly or extracted from the analysis of i3 — > D^*^D*q{Dsi) decays in 
Refs. Pfil'&E [Tl |. [l^l-fl^. The form factors of _B ^ D^*^lvi transitions have been analyzed and calculated in 
various theoretical approaches such as: heavy quark effective theory, QCD sum rules, lattice QCD, different types of 
quark and soliton models, approaches based on the solution of Bethe-Salpeter and Faddeev equations, etc. 

In this paper we assume the D*^ and Dgi mesons to be hadronic molecules - bound states of D,K and D*,K 
mesons, respectively. Using this molecular picture, in Refs. [H, |23| we calculated strong and radiative decays of D^q 
and Dsi mesons. The obtained results are in agreement with other theoretical approaches, e.g. the strong decay 
widths D*Q DsTT^ and Dsi -D*7r° are of the order of 10^ KeV and the radiative decays D*q D*j, Dgi Dgj, 
etc. are of the order of a few KeV. Here, using the same approach, we calculate the leptonic decay constants /d*^ and 
fosi- For the form factors governing the semileptonic B D'-*'>£Di transitions we use the model-independent results 
obtained by Caprini, Lellouch and Neubert (CLN) ''2^ on the basis of heavy quark spin symmetry and dispersive 
constraints, including short-distance and power corrections. Note, that in Ref. ^18,] the authors already used the CLN 
results in their analysis of two-body B DD*q{Dsi) transitions restricting themselves to the heavy quark limit and 
modes with a pseudoscalar D meson in the final state. Using the experimental lower limits for the B — DD*q{Dsi) 
branching ratios they derived lower limits for the products |ai|/£)»^ and |ai|/Dsi, where oi is a combination of the 
short-distance Wilson coefficients (29j-[3]|. 

In the present paper we proceed as follows. First, in Section II, we discuss the basic notions of our approach. 
We indicate and evaluate the effective mesonic Lagrangian for the treatment of charmed mesons I?*o(2317) and 
Dsi (2460) as DK and D*K bound states, respectively. Then in Section III we discuss the calculation of the leptonic 
decay constants /d*^ and fosi- In Section IV we present a detailed analysis of two-body bottom meson decays 
B — + D'^*^ DIq{Dsi) applying the factorization hypothesis. As we already stressed before, in this analysis we use the 
model-independent CLN results [2^ for the weak form factors defining the B — > D^*^ivi transitions. In Section V we 
give a short summary of our results. 



II. MOLECULAR STRUCTURE OF ^^(±(2317) AND D±(2460) MESONS 

In this section we discuss the formalism for the study of the 15*0^(2317) and D^^{2A6Q) mesons as hadronic molecules, 
represented by DK and D*K bound states, respectively. We adopt that the isospin, spin and parity quantum numbers 
of D;o^(2317) and Dfi(2460) are: /(J^) = 0(0+) and /(J^) = 0(1+), while for their masses we take the values: 
TO£).^ = 2.3173 GeV and m£)^^ = 2.4589 GeV |1(|. Our framework is based on effective interaction Lagrangians 
describing the couphngs of 15*0(2317) and Dsi(2460) mesons to their constituents: 

^D-'Jx) =5^.^d:o-(x) Jdyl>D'Jy^)D{x + w^^y)K{x~w^^y) + B.C., (1) 

^dA^) = 9n^,D^-{x) fdy^DAy')D;{x + w,,,,y)K{x-w,,,y) +Ii.c., (2) 
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where the doublets of D^*^ and K mesons are defined as 



D = 



D+ 



D* = 



D* + 



K = 



K+ 



(3) 



The summation over isospin indices is understood. The molecular structure of the -D*^ and Df^ states is: 



1+^ 



\d+k") 



\D:,) = \D*+K') 



\d"k+) , 
\d*"k+) 



\D 



sO 



= \D-K 



0\ 



|i?7i) - \D*-K 



0\ 



(4) 



The correlation functions $m with M — D*q or Dsi characterize the finite size of the Z?*q(2317) and Dsi{2A60) 
mesons as DK and D*K bound states and depend on the relative Jacobi coordinate y with, in addition, x being 
the center of mass (CM) coordinate. Note, that the local limit corresponds to the substitution of by the Dirac 
delta-function: — > 5'^{y)- In Eqs. ^ and ([2]) we introduced the kinematical parameters w^ : 



m,i 



Wi. 



rrii + m,- 



(5) 



where mo, "niD- and ttik are the masses of D, D* and K mesons. The Fourier transform of the correlation function 
reads 



$m(V), M^D:^,Dsi. 



(6) 



A basic requirement for the choice of an explicit form of the correlation function is that it falls down sufficiently fast 
in the ultraviolet region of Euclidean space to render the Feynman diagrams ultraviolet finite. We adopt the Gaussian 
form 



^m{p%) = exp(-p|/A 



M) 7 



(7) 



for the vertex function, where pE is the Euclidean Jacobi momentum. Here A/^*^ is a size parameter, which 
parametrizes the distribution of D and K mesons inside the D*^ molecule, while A^^j is the size parameter for 
the Dal molecule. For simplicity we will use a universal scale parameter Am = Ao*^ = A^^j^, i.e. the same size for 
D*o and D^i mesons. 

The coupling constants 50*^ and go^i are determined by the compositeness condition (3^. [33j. which implies that 
the renormalization constant of the hadron wave function is set equal to zero: 



1-Ei5.^(ml,.j = 0, 



(8) 
(9) 



Here, S'^. (™l3» ) = 3^, n^. {m\,, ) is the derivative of the D*^ meson mass operator. In the case of the Dgi meson 
we have 1) = ffo II'd 1)1 which is the derivative of the transverse part of its mass operator S^''^, 

conventionally split into transverse "^d^i and longitudinal 5^13^^ parts as: 



,2 ^D, 



(10) 



where 



91 



9' 



p^'p" 



ay 

9± 



0. 



The mass operators of the D*g and Dsi mesons are described by the diagram of Fig.l. 

Following Eqs. ((8]) and ^ the coupling constants g^. and ^ can be expressed in the form: 



(11) 



1 



(47rAM)2 
2 

(47rAM)2 



1 00 
dx 







dx 



da a Po{a, x) 
(l + a)3 

daaPi{a,x) 
(l + a)3 



2^|j.(l + a) dzi 



*m(^i) 



(12) 
(13) 
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where 



Po{a,x) = a'^x{l — x) + w^^^ax + w'^^^a{l ~ x) , 

Pi{a,x) — a'^x{l — x) + w'^ ^ ax + w'^ ,a(l — x), (14) 



zo = ^^aa; + ^^a(l -x) — fx^. 

1 + a "0 

2 ,2/1 N -Pi (0^,2:^) 2 

I + a 



^J■M 



Above expressions are valid for any functional form of the correlation function <i>M . 

Let us note, that the compositeness condition of the type (H]), ([9]) was originally applied to the study of the 
deuteron as a bound state of proton and neutron (3^ . Then this condition was extensively used in low-energy hadron 
phenomenology as the master equation for the treatment of mesons and baryons as bound states of light and heavy 
constituent quarks |33|-[^. In Refs. [26|, |3^, [s^l this condition was used in the application to glueballs as bound 
states of gluons and light and heavy scalar mesons as hadronic molecules. 



III. LEPTONIC DECAY CONSTANTS AND /„ 

The decay constants of the scalar £)*q and axial Dsi mesons are defined by 

{Q\-s0^c\D:+{p))=p^f^,^, (15) 
(0|sO^c|i?+(p,e)) = e^{p)mu,J,^^ , (16) 

where — — 75). Note, parity symmetry implies that only the vector component of the weak V — A current 
contributes to the transition D*^ W~^, while for the transition D'^^ only the axial component is present. 

The one-loop meson diagrams describing the matrix elements of Eqs. p5|) . ^6]) are given in Figs. 2a and 2b. In other 
words, adopting the molecular picture for the D*q and D^i mesons we need to evaluate the two-point DK loop diagram 
describing the transition of the D*q meson to the vector current and the two-point D*K loop diagram corresponds to 
the transition of the D^i meson to the axial current. The coupling of the weak c — s- s flavor-changing vector current 
(FCVC) Vfj. to the DK pair and the weak c — > s flavor-changing axial current (FCAC) Afj_ to the D*K pair can be 
extracted from data. In particular, the couplings are extracted by deriving the effective Lagrangians from the matrix 
elements for the D — > K^*^lvi semileptonic transitions. First, the matrix element for the D K£vi transition ^sl] 
at zero momentum transfer g = p — p' — )■ is approximated by 

{K{p')\v^mD{p)) f^'H^){p+pr ■ (17) 

Therefore, the effective Lagrangian describing the coupling of D and K mesons with the weak vector current is given 
by 

Cvdk{x) = /^^(O) y'^(x) D{x) td- K{x) -f H.c. (18) 

where Ad-B = (9^/1)5 ~ Adf,B. 

The effective coupling of the weak axial current A^ with D* and K mesons can be related to the corresponding 
coupling with D and K* mesons on the basis of SU(4) flavor symmetry arguments. Both matrix elements for the 
semileptonic transitions D* KiD^ and D —f K*£i'£ semileptonic transitions at zero momentum transfer can be 
approximately written as: 

{K{p')\A^{0)\D*{p,e)) 6^'(mz,.+mA')Af ^(0), (19) 
{K*{p',e')\A'^mD{p)) ^ e'*^imD+fnK')Af''\0). (20) 

At present no experimental information on the axial form factor A^ ^(0) is available, while the form factor A^^ (0) 
is known, although with significant uncertainties. SU(4) symmetry requires that the axial form factors ^(0) and 
A'^'^' (0) satisfy the relation: 



imD' + mK)Af'^0) = (niD + mK')A^'^' (0) , 



(21) 
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which in turn allows to express the unknown form factor ^ (0) by Af^ (0). The effective Lagrangian, describing 
the coupling of D*K and DK* meson pairs to the weak c ^ s FCAC is then given as: 



Cad'k{x) + Cadk' (x) = {mo + mK- ) A^"^' (0) A^(x)\dI{x)K{x) + D(x)Kl{x) 



H.. 



In the following we calculate the leptonic decay constants / and /„ based on the effective interaction Lagrangian 

£eff, which includes the couplings of D*o, -Dsi mesons and the weak currents (vector and axial A^) with DK and 
D*K meson pairs: 



(23) 



The corresponding diagrams are given in Fig. 2. After a straightforward evaluation we obtain the following analytical 
expressions for the leptonic decay constants: 



9o' 



8tt 



da a 
[l + a)' 



2{w^^ + ax) 



(24) 



87r2 



mo, I 



A\"^ (0) / dx 



da a 
(l + a)- 



a dzi 
ifijj, da 



(25) 



where the coupling constants g^, and g^^^ are defined by Eqs. (fT2|) and (fT3|) . We have only one model parameter 

in our calculations - the scale parameter Am which was previously fixed to the value of Am = 2 GeV [1^ [13] from 
strong and radiative decays of D*q and D^i mesons. Further quantities, entering in Eqs. ([M)) and (P5|) . are chosen as 
follows. The masses of mesons we take from the PDG USM 



m_D*^ = 2.3173 GeV , "t-d^i 
rriK = mK± = 493.677 MeV , m^. 
r7i£i = mjji — 1.8693 GeV , too. 



= 2.4589 GeV , 
: mK,± = 891.66 MeV, 
TOc-i = 2.010 GeV. 



(26) 



For the values of the D - 
world average data [H, Ho 



K^*'>£i'e semileptonic form factors at zero recoil, entering in our calculation, we use the 
m of: 



f'^^iO) = 0.75 ± 0.05 , ^* (0) = 0.65 ± 0.05 . 

Using this input we get the following results for f^, and f^^ ^ : 

= 67.1 ±4.5 MeV, /^^^ = 144.5 ± 11.1 MeV . 



(27) 



(28) 



In Table 1 we summarize the present results for f^, and /j-, ^ obtained in different approaches (either on the basis 

of hadronic models or from the analysis of exp erimental data on two-body i3-meson decays). Our results are in 
agreement with the predictions of Refs. p^.ll7l.[lS, 2C|, especially with the lower limits derived from an analysis of 
the branching ratios of B ^ D^*^DIq{Dsi) decayslllTSl . 



At this point it is worth discussing the heavy quark limit (HQL) for the leptonic decay constants f^, and /j-, ^ , 
where the masses of D, D*, D*q and Dsi mesons together with the charm quark mass TOc approach infinity. In the 
HQL the D'^*^ mesons in the D*q{Dsi) hadronic molecules move to the center of mass and are surrounded by a light 
K meson in analogy with the heavy-light Qq mesons. It is known (see e.g. discussion in Refs. [l7|, Esl) that in the 
two-quark picture for the D*^, Dgi mesons HQL leads to degenerate values for these couplings 

fo. ^ fo^, - ■ (29) 



In the present molecular approach we can also guarantee that the couplings are degenerate in the HQL, but the scaling 
law is different from the behavior. In a first step we apply the HQL to the coupling constants g^, and g^^ ^ , 

which are degenerate: 

1 _ 1 lo 
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/n = 



5 



da 



(47rTOc)^ 



(30) 
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1 OO 

The structure integrals J dx J da - • ■ entering in Eqs. ([^ and (^51) are also degenerate and equal 



OO 

Am f da y/a 2 f 2 /oia 
rric J 1 + a 



Finally, the leptonic decay constants and ^ in the HQL are given by: 

_ hql J±_ /on) 



where we introduced the effective VDK and AD*K couphngs in the HQL: g'^^'^j^^ and ff^^'.^^- From Eqs. ([5^ and 
we deduce that in the HQL / and /„ in the HQL limit do not depend on rric at all, unlike as in Eq. (I^S)) . and 

are degenerate for g^!;'^'^ = ff^^'.^, • The latter condition can be eventually fulfilled, e.g. at finite masses (see Eqs. ([M)) 
and (|25|l) the ratio gAo^fc/dvoK is close to 1: 

To give an estimate for the absolute values of the leptonic decay constants f^^, and ^ in the HQL we vary the 
coupling constants g^''^^ = g^j^^^^^ in the region 0.75 ± 0.25. The corresponding result is: 

= /^^^ = 205.2 ± 68.4 MeV . (35) 

The following comment is in order. As was already stressed in Rcf. ,18.j the HQL does not give a reasonable approx- 
imation for the P-wave Dgj meson system. According to the data on the two-body decays B D^*^DIq{Dsi) the 
physical value of the decay constant ^ should be about twice as large as f^, [isl. [20l|. 



IV. WEAK TWO-BODY DECAYS B D(*'D*o(-Dsi) 

In this section we give the predictions for the branching ratios of i? ^ D^*^ D*q{Dsi) decays. For this purpose we use 
leptonic decay constants , ^ and, in addition, model- independent results for the form factors of i? ^ D^*^ivi 

transitions obtained by Caprini, Lellouch and Neubert (CLN) [ll]. Latter derivations are based on heav y q uark 
spin symmetry, dispersive constraints, including short-distance and power corrections. Note, that in Ref. |l8| the 
authors already used the CLN results in their analysis of two-body B DD*q{Dsi) transitions, restricting to modes 
containing the pseudoscalar meson D in the final state. 

Working with the factorization approximation we first write down the factorizable amplitudes for the four decay 
modes D^D*^ , D^D^^, D*^D*q and D*^D^^ as expressed in terms of matrix elements of the semileptonic 

B — > D^*Hvi and the leptonic D^q, Dgi transitions and (see details in Refs. pTl [Tsl. [20| ) : 

M{B^^D+d:^) = G,n{D:„-{q)\rsO^c\0) {D+{p')\cO'^b\B%p)) , (36) 

M{B°^D*+d:^) = G,s{D:,-{q)\sO^c\0) {D*+{p',e')\m\B''{p)) , (37) 

M{B°^D+D;,) = G,s{D;,{q,e)\sO^c\0) {D+{p)\cO'^b\B°{p)) , (38) 

M{B"^D*+D;,) = G,s{D;,{q,^)\sO^c\0){D*+{p',e')\m\B"{p)). (39) 

Here 

Geff = ^K6K>1 (40) 

and fli = C2 + ci/Nc is the combination of the short-distance Wilson coefficients ci and C2 [2^-[3l|. Using the 
Wirbel-Stech-Bauer (WSB) decomposition [Hjli^l of the hadronic B — > Z?'*) matrix elements 

{D+ip')\V^\B\p)) = (^{p + pT ~ ""^ ",""^ ^"1^1(9^)+ ""^ ~,"''^ g"fo(g^), (41) 
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{D*+{p',e')\V^\B\p)) = e':p'^p,V{q^), (42) 

ruB + ruD' 

{D*+{p',t')\A^^\B\p)) = e'*^{mB+mD^)A,{q^) 

{p + prMq') + 2rnn^q^'-^Uo{q')-As{q')}, (43) 



where — c^'^b, — c^'^^^b and 

Ml ) = — 7y ^iW ) 7y Mq ) (44) 

we arrive at [13, [H, [ll : 

M{B° ^ D+D;-) = Geff /^.^ (m^j - ml) F^{ml,J , (45) 

MiB'^^D+D-^) = G,sfn^,mD^,e*ip + p')Fi{mlJ, (46) 

M(S°^i?*+i?:o-) = 2G,sf„,^mD'e'*pAo{ml.J, (47) 

M{B° - = Geff /^^^ mz,,, e;e':p'^pf, VimlJ 

[mB+mD' 

+ e'*e*{mB + mD*)A,{mlJ- ^''*P'*P A^imlj] . (48) 

TUB + niD' J 

In the following it is convenient to express the WSB set of form factors through a set of form factors hjiw) depending 
on the kinematical variable w — v ■ v' , the scalar product of the four-velocities of B and D'-*^ mesons 0, |42|. [isj: 

{D{v')\V^'\B{v)) = h+{w){v + v'Y + h_{w){v-v'Y , (49) 
{D*{v',e')\V>^\B{v)) = ihv{w)€>^''"^t':v'^vp, (50) 



{D*{v\e')\A^^\B{v)) = hA^w) {w + l)e'*^ 



(51) 



where the meson states \M{v)) obey the mass-independent normalization condition 

{M{v')\M{v)) = 2t."(27r)353(p _ p') ^ (52) 

instead of the relativistic one used before: 

{M{p')\M{p)) - 2/(27r)353(p _ p') . (53) 

The WSB form factors are expressed in terms of the hi{w) form factors as: 



F^{q^) = ^V^iw), (54) 



2, 2^ w + 1 



Fo{q')^j^^V,{w)R{w), (55) 



2^ 1 + ^* 



Viq') = —^hAAw)Riiw), (56) 
2\/r* 



2, 2^/r* w + 1 



A^iq') = ^^^-hAAw), (57) 
A2{q^) = hA,H R2{w) , (58) 



2^^* 
1 

2Vr* ' 



Ao{q^)^ -^hAAw)RoH, (59) 



where r = mu/mB, r* = mD*/mB and 

l + r 



Vi{w)=G{w)^h+{w)-\ — -h_{w). (60) 
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Here we use the well-known form factor ratios i?i (w) and R2 (w) [2^, ]^ : 

hyjw) hA^jw) +r*hA2{w) 
-^iW^I — TT' R2(w) = — ^ . 61) 

In the literature the form factor Vi{w) is also denoted as G{w). In order to express all hi{w) form factors, completely 
defining the B D^*^ transitions, through the two form factors Vi{w) and hAi{w) we introduce the additional ratios 
R{w), R*2{w): 

1 -\- r w — 1 

h+iw) -h-(w) f, / ^ *j, / ^ 

u/ \ 1-rw + l ^ ' r}*( ^ hA^{w) -r*hAAw) , . 

R{w) = r—p , ^2 H = 7 — r~\ • 



and Rq{w), which is just the combination of R2{w) and R2{w): 



Roiw) = {w + l)[l-R;{w)] + 



(l + r*)2 



2r* 



R;iw)-^-^R2{w) 
1 + r* 



(63) 



For the functions Vi{w) = G{w), hAiiw), Ro{w), Ri{w), i?2(w) and i?2(^) ^^e the model-independent results 
derived by Caprini, Lellouch and Neubert (CLN) [28]: 

^^'^^ - 1 - 8p%z + {51p% - lQ)z^ - (252/9^ - 84)^^ , (64) 



G(l) 

'"^ 1 - z + (^ipl - I5)z2 _ (231p2 _ 91)^3 ^ (g5) 



i?(w) = 1.004 - 0.007(iy - 1) + 0.002(w - 1)^ , (66) 

i?i(w) = 1.27 - 0.12(w - 1) + 0.05(w - 1)^ , (67) 

i?2(w) = 0.80 + 0.11(?i; - 1) - 0.06(w - 1)^ , (68) 

K*2{w) = 1.15 - 0.07(w - 1) - 0.11(?i; - 1)^ , (69) 



where 



z = — . (70) 

In the derivation of the expressions for the ratios Rq {w) and i?2 {w) we used the results of Ref. [2^ (see Appendix) . 
Note, that R{w) is very close to 1 for w in the interval 1 < w < lUmax = + Tn'jj)/{2mBmu). 

The two-body decay widths T{B — > D^*^D*g) and T{B D'^*^D^^), given in terms of the CLN form factors, are 
expressed by the formulas (see e.g. Ref. 18]): 



r(i? - 




T{B - 


- DD^i) = 


T{B - 




r{B - 





— (ms - mof ml {wi + if J wj - 1 [G{wi) R{wi)f (71) 
+ mDfml{wl-lf^^G^{w2) (72) 



- ^ fl, < i^l - ^f'^ [hA. M RoM? (73) 



where the kinematical variables w, are defined as follows 



- mo' )^ ml. {wi + If ^f4~l h\ (wi) Pha, (wi) , (74) 



ml + ml -ml. ml + ml - ml^ 



= 2^^B ' 2mBmn ' 

ml + ml. - ml. ml + ml. - m^, 

W3 = ^ , W4 = , (75) 

ZmBmD' ImBmn* 
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and 



1 — 2wr* + r 
' (l-r*)2 



*2 



W — 1 
W+1 



W — 1 

1 - r* 



(l-i?2M) 



(76) 



Note that the product h\_^('w) Pha^ (w) defines the well-known function F{w) [H, E^l, which governs the semileptonic 
transition B — > _D*i?p£ (see e.g. Eq.(35) in Ref. psj): 



1 



Aw 1 - 2wr* + r* 



w + 1 (l-r*)2 



(77) 



In the numerical calculation we use the value for the Cabibbo-Kobayashi-Maskawa (CKM) matrix element Vcs = 
0.97296 ± 0.00024 from the 2006 global fit [10, Sll and the 2006 averaged values from the Heavy Flavor Averaging 
Group (HFAG) [ii]: 



\Vcb\F{l) = (36.2 ±0.8) X IQ-^ 
\Vcb\G{l) = (42.4 ±4.4) X 10"^ 



pI = 1.19 ±0.06, 
p% = 1.17±0.18, 



(78) 
(79) 



where the normalization of F{w) and its slope a,t w — 1 are related to the characteristics of the /iai('w) form factor 
as 



Ph, 



0.21. 



(80) 



For the parameter oi we use the value of 1.05 from Ref. [3l|. A detailed discussion concerning the choice of ai can 
be found in Ref. [30| . 

In Table 2 we present our predictions for the branching ratios of two-body decays B D^*^DIq{Dsi). For the data 
we use the averaged lower limits from PDG [l3] and in addition for the modes with £'si(2460) meson in the final state 
the direct results of the BABAR Collaboration Our predictions are in good agreement with the experimental 
data except for the marginal situation in the case of B° D*q D*~^ decay, where our prediction is slightly lower 
than the experimental limit. In Table 3 we present the results for the ratios r{B D* D*Q)/r{B — > DD*q) and 
T{B D* Dsi) /T{B DDsi). We also use the compilation of experimental data and theoretical results within 
the covariant light-front (CLF) approach [l3| summarized in Table 10 of Ref. Here, our predictions are in good 
agreement with the existing experimental data. In comparison with the CLF approach our predictions are lower, 
although not significantly. 



V. CONCLUSION 



We considered the new charmed-strange mesons 15*0(2317) and 1551(2460) as hadronic molecules, that is DK 
and D*K bound states, respectively. Using an effective Lagrangian approach describing the coupling of _D*q(2317) 
and _Dsi(2460) states with their constituents we determined the leptonic decay constants fo'^ and fo^i- Then we 
presented a detailed analysis of two-body bottom mesons decays B ^ D*^*^ D;„{D si) using the factorization hypothesis 
and model-independent results for the form factors of the B D^*HPe transitions obtained by Caprini, Lellouch and 
Neubert (CLN) [2^. The decay widths are derived in terms of the CLN form factors. Calculated branching ratios 
for B D^*'^D*q{D si) decays are in good agreement with experimental results (or the lower limits), supporting a 
possible interpretation of D*q(2317) and Dsi{2A60) as hadronic molecules. 
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Table 1. Leptonic decay constants f^, and /j 



Approach 


L> (MeV) 


fo,, (MeV) 


Rcf. [23] 


225 ± 25 


225 ± 25 


Ref. [24] 


206 ± 120 




Rcf. [25] 


200±50 




Ref. [21] 


170 ± 20 


247 ± 37 


Ref. [14] 


138±16 


259±13 


Ref. [22] 


110±18 


233±31 


Ref. [18] 


> (74 ± ll)/|ai| 


> (166 ± 20)/lail 


Ref. [17] 


71 


117 


Ref. [17] 


60± 13 


150 ± 40 


Ref. [20] 


> (58 - 86)/lail 


> (90 - 228)/lail 


Ref. [12] 


67±13 




Ref. [11] 


44 


41 


Our results 


67.1 ± 4.5 


144.5 ± 11.1 



Table 2. Branching ratios of _B ^ D'-*^ DIq{Dsi) decays (in units of 10 ^). 



Mode 


Data (averaged) [10] 


BABAR [9] 


Our results 


B- ^ d:^d" 


> n 74+0-23 




1.03 ± 0.14 




> 0.97±°1i 




0.96 ± 0.13 


B- ^ d;,d" 


> 1.4t^J 


4.3 ± 1.6 ± 1.3 


2.54 ±0.39 


BO ^ d;,d+ 


> 


2.6 ±1.5 ±0.74 


2.36 ± 0.36 


B- ^ d:oD*o 


> 0.9 ± o.eto.t 




0.50 ± 0.07 


B° D*-D*+ 


> 1.5±0.4l°j 




0.47 ± 0.06 


B- D^^D*" 


> 5.5± 1.2l? j7.6 


11.2±2.6±2.0 


7.33 ± 1.12 


B° D-^D*+ 


> 7.6± 1.7l|4 


8.8 ± 2.0 ± 1.4 


6.85 ± 1.05 



Table 3. Ratios MD*/MD = T{B ^ MD*)/T{B MD) for M ^ D*q, D 



Mode 


Data [20] 


CLF [20] 


Our results 




0.91 ±0.73 


0.49 


0.48 


d:^d*+/d:,d+ 


0.59 ±0.26 


0.49 


0.48 




3.4 ±2.4 


3.6 


2.9 


d;,d*+/d;,d+ 


2.6 ± 1.5 


3.6 


2.9 
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FIG. 2: Diagrams related to the leptonic decay constants of Z)*o(2317) and Dsi(2460) mesons. 



